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We develop the real stabilization method within the framework of the relativistic mean field 
(RMF) model. With the self-consistent nuclear potentials from the RMF model, the real stabilization 
method is used to study single-particle resonant states in spherical nuclei. As examples, the energies, 
widths and wave functions of low-lying neutron resonant states in !?°Sn are obtained. These results 
are compared with those from the scattering phase shift method and the analytic continuation in 
the coupling constant approach and satisfactory agreements are found. 
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I. INTRODUCTION 


The investigation of continuum and resonant states is an important subject in quantum physics. In recent years, 
there has been an increasing interest in the exploration of nuclear single particle states in the continuum. The 
construction of the radioactive ion beam facilities makes it possible to study exotic nuclei with unusual N/Z ratios. 
In these nuclei, the Fermi surface is usually close to the particle continuum, thus the contribution of the continuum 
and/or resonances being essential for exotic nuclear phenomena [1, 2, 3, 4, 5, 6]. It has been also revealed that the 
contribution of the continuum to the giant resonances mainly comes from single-particle resonant states [7, 8]. 

For the theoretical determination of resonant parameters (the energy and the width), several bound-state-like 
methods have been developed. The complex scaling method (CSM) describes the discrete bound and resonant states 
on the same footing [9]. In this method, a complex coordinate scaling is introduced to rotate the continuum into 
the complex energy plane and the wave functions of resonant states, but not scattering states, are transformed into 
square-integrable functions [10]. Although it involves the solution of a complex eigenvalue problem which causes some 
difficulties in practice, the CSM has been widely and successfully used to study resonances in atomic and molecular 
systems |11, 12, 13] and atomic nuclei [9, 10, 14, 15, 16]. The analytical continuation in the coupling constant (ACCC) 
approach is based on an intuitive idea that a resonant state can be lowered to be bound when the potential becomes 
more attractive or equivalently the coupling constant stronger, thus a resonant state being related to a series of 
bound states via an analytical continuation in the coupling constant [17, 18, 19]. Combined with the cluster model, 
the ACCC approach has been used to calculate the resonant energies and widths in some light nuclei [20, 21]. An 
attempt to explore the unbound states by the ACCC approach within relativistic mean field (RMF) model was first 
made in Ref. [22] where resonant parameters of some low lying resonant states in 160 and “°Ca obtained from the 
ACCC calculations are comparable with available data. The wave functions of nuclear resonant states were also 
determined by the ACCC method where the bound states are obtained by solving either the Schrödinger equation 
with a Woods-Saxon potential [23] or the Dirac equation with self-consistent RMF potentials [24]. 

The real stabilization method (RSM) is another bound-state-like method [25]. The equation of motion of the system 
in question is solved in a basis [25] or a box [26] of finite sizes, thus a bound state problem being always imposed. The 
RSM uses the fact that the energy of a “resonant” state is stable against changes of the sizes of the basis or the box. 
It has been used to calculate the resonance parameters in elastic and inelastic scattering processes [27, 28]. Some 
efforts have also been made in order to calculate more efficiently resonance parameters with the RSM [29, 30, 31, 32]. 
In this work, we investigate single particle resonances in atomic nuclei by combining the RSM and the relativistic 
mean field (RMF) model. 

The paper is organized as follows. In Sec. II we give briefly the formalism for the RSM and the RMF model. The 
numerical details, the results for !2°Sn and discussions are given in Sec. III. Finally we give a brief summary. 
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II. FORMALISM OF THE RMF MODEL AND THE RSM 
A. The relativistic mean field model 


The basic ansatz of the relativistic mean field (RMF) model is a Lagrangian density where nucleons are described 
as Dirac spinors which interact via the exchange of several mesons (a, w, and p) and the photon [33, 34, 35, 36, 37], 


L = Filip- M) i + 59,00"0 — U(0) — gofio 
Lace + iniu — gu 
= FRB” + S270" — oohi BTh 
= TF Fe = eji Adi, (1) 
where the summation convention is used and the summation over 7 runs over all nucleons, ¢ = y#x, = yx", M the 


nucleon mass, and Mo, Jo, Mw, Jw, Mp, Jp Masses and coupling constants of the respective mesons. The nonlinear 
self-coupling for the scalar mesons is given by [38] 


U(o) = =m2o07 + P4 Ba, (2) 


and field tensors for the vector mesons and the photon fields are defined as 


iw = OMEP ~ OW, 
Ruy = ppo — WP — Gp(Pu X Pv), (3) 
Fup = OY T Op Ap. 


The classical variation principle gives equations of motion for the nucleon, mesons and the photon. As in many 
applications, we study the ground state properties of nuclei with time reversal symmetry, thus the nucleon spinors 
are the eigenvectors of the stationary Dirac equation 


[a-p+V(r) +M + S(r))] dir) = cedi(r), (4) 


and equations of motion for mesons and the photon are 


n E a = a 

= +m? w (r) = gupr(r), 

(-A +m?) p°(r) = gppa(r), (5) 
—AA(r) = ep(r), 


where w? and A? are time-like components of the vector w and the photon fields and p° the 3-component of the 
time-like component of the iso-vector vector meson. Equations (4) and (5) are coupled to each other by the vector 
and scalar potentials 


1 
V(r) = gu (r) + gprs? (r) +e A.r), © 
S(r) = goo(r), 
and various densities 
psr) = Di Gilt Wil), 
plr) = diel wl (r)bi(r), 7 
plr) = Ei V(b), 5 
—T: 
pel) = Dies Bi (r) tlr). 
For spherical nuclei, meson fields and densities depend only on the radial coordinate r, the Dirac spinor reads 
GE(r 
Ce yr 4, 4, s) . 1 
Parm(r, 8, t) = k Xta (t), J = It 5? (8) 


FE(r) i 
EO yt. Q, s) 
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with Yhn (0, p) the spin spherical harmonics. The radial equation of the Dirac spinor, Eq. (4), is reduced as 


as = (-2 4%) FEAM HSV, 

Or r (9) 
eben EE E E 

ðr r 


The meson field equations become simply radial Laplace equations of the form 


o> 20 2 
-Z -rar m3) olr) = sg(r). (10) 
mg are the meson masses for ¢ = g,w, p and zero for the photon. The source terms are 


—JoPs(r) — g20°(r) = g30?(r), for o, 


z Gapa), for W, 
900) = ) goal), for p, (1) 
epe(r), for A, 
with 
4nr’p.(r) = Die (IGi(r)? -= |Fi(r)|?), 
4nr*py(r) = Lier (IGi(r)/? + |Fi(r)|?), (12) 
Arr? p3(r) = X; 2ti(|Gi(r) |? + |Fi(r)|7), 
Anr*pe(r) = Yia (3 — ti) (IGi(r)? + |F). 


The above coupled equations can be solved iteratively in r space [39] or in the harmonic oscillator basis [40] using the 
no sea and the mean field approximations. 


B. The real stabilization method in coordinate space 


With the self consistent vector and scalar potentials V(r) and S(r), the Dirac equation (9) is solved in a spherical 
box of the size Rmax under the box boundary condition, and thus the continuum is discretized. When Rmax is large 
enough, the energy of a bound state does not change with Rmax. In the continuum region, there are some states 
stable against the size of the box, i.e., the energy of each of such states is almost constant with changing Rmax; such 
stable states correspond to resonances. 

The resonant parameters, Æ, and I’, may be obtained by fitting the energy E and the phase shift 7(£) in an energy 
range around a resonance to the following formula [25], 


m(E) = npor(E) + tan (5) (13) 


The phase shift m(E) can be calculated as [41] 


Rmax ; 
n ( z) oo fi xi(r) [E — H(r)] f(r) sin krdr (14) 


Le ’ 

"o3 fo" alr) [E — H(r)] f(r) cos krdr 
with f(r) satisfying f(r) — 1 when r > oo and f(0) = f’(0) = 0. However, Eq. (14) converges very slowly with the 
box size due to the influence of the non zero centrifugal potential at large Rmax [41]. 

In the present work, we use a simpler method proposed by Maier et al. [26] in which it’s not necessary to calculate 
the phase shift. The resonance energy is determined by the condition 6?E/3R? ax = 0 and the corresponding box 
size is labeled as Rmax, i.€., Ey,=E (Rmax). The width is evaluated from the stability behavior of the positive energy 
state against the box size around Rmax- 

When r is large enough, the nuclear potentials S(r) and V(r) vanish, G(r)/r satisfy 


HG | (e — a G=0, i) 


dr? r2 
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with a? = E? — M?. The general solution reads 
G(r) x ar [cosm ji(ar) — sin m ni(ar)] , (16) 
When r — œ, G(r) x sin(ar — m +m). Therefore when Rmax is large enough, 


l 
aRmax = = +H = nT. (17) 
Under the assumption that the phase shift from the potential scattering m,p,o1(E) varies slowly with respect to the 
box size, ie, Om pot/ORmax ~ 0, one derives from Eqs. (13) and (17) the formula, 


F 2,/ E2 + 2EM 


A (18) 
—(Ey + M)Rmax — (E2 + 2EM) [9E/ORmaxlk,,..] 


In the non-relativistic limit, Ky « M, Eq. (18) is reduced to 


"N 2,/2E,/M (19) 


—Rinax = 2Ey [dE /dRmaxl Rna] = 


which is essentially the same as Eq. (8) in Ref. [26] except for that here natural units with A = c = 1 is used. 


II. RESULTS AND DISSCUSION 


In this section we present the results of the RSM in the framework of the RMF model. In our calculation, we use 
for the Lagrangian density the effective interactions PK1 [42] and NL3 [43]. We change the size of the box in a large 
range (7 fm < Rmax < 60 fm if not specified) in order to find not only narrow resonances but also wide ones. We 
take 1?°Sn as an example and compare the results for neutron resonances from the RSM with those from the ACCC 
approach [24] and the scattering phase shift method [44]. 

By examining the stability of low lying positive energy states against the size of the box Rmax, we find that except 
for s state, there are neutron resonances in !?°Sn with the orbital angular momentum l up to 6. 

The narrowest resonance is an i,3/2 state which lies at about 3.45 MeV above the threshold. The positive energy 
Vi13/2 States in a box of different sizes are shown in Fig. 1. With the box size Rmax increasing, the lowest vi;3/2 state 
first falls down quickly then its energy becomes constant in a large region of Rmax. After it crosses with the second 
lowest viis /2 State at around Rmax = 25 fm, the lowest one falls down again. Similar level-crossings occur regularly 
at larger Rmax. This stability behavior implies that there is a narrow vi,3/2 resonant state with approximate energy 
3.45 MeV. The resonant energy 3.469 MeV and Rmax = 20.1 fm are obtained under the condition 0?E/OR?,,,. = 0 


max 


neutron 


FIG. 1: Positive energy vij3/2 states in 20Sn under different box boundary conditions. 


201708.00263v1 


chinaXiv 


0.2 


fo) 
2x 


S 
Z 


E =3.469MeV 
T =0.005MeV 


. -3/2 
wave function(fm™ į) 
=] 
oO 


ò 
N 


FIG. 2: The radial wave function for the neutron vij3/2 resonant state in 1205p, 


TABLE I: Energies and widths of the single neutron resonant state i13/2 in 1205) from different Æ ~ Rmax curves (cf. Fig. 1). 
Rmax is in fm and E} and T are in MeV. 
Rmax 


E POTR] e OT 
20.1 |3.4688 | 0.327 || 47.2 |3.4686 | 0.470 
30.3 13.4686 | 0.426 || 55.3 13.4685 | 0.478 
38.9 [3-4686] 0.456 |] 63.2 [3.4686 | 0.483_ 9 13.4686] 0.456 || 63.2 |3.4686] 0.483 


from the first E ~ Rmax curve (labeled as “1i,3/2”) in Fig. 1. The width 0.003 MeV are obtained from Eq. (18). One 
of the approximations made in deriving Eq. (18) is Rmax should be large. We next examine the dependence of the 
resonant parameters on the box size by calculating FE, and [ from other E ~ Rmax curves with larger Rmax. For 
this purpose the calculations with Rmax up to 65 fm are carried out. The results are given in Table I. The energy is 
almost a constant with increasing Rmax. The variation between the widths obtained from adjacent E ~ Rmax Curves 
decreases with Rmax and is about 1% at Rmax ~ 63 fm. For other resonant states presented in this work, we also 
make similar investigations. Once it converges to within 1%, the value of the width is assigned to a resonant state. 
The wave function of the resonant state vi13/2 is given in Fig. 2. From this figure one can find that this state is almost 
localized inside the nucleus which is consistent with the small width rI. 

Although it lies below the state vi13/2, the resonant state vfs/2 is about an order of magnitude wider than vi15/2- 


FIG. 3: Positive energy vfs/2 states in 120S) under different box boundary conditions. 
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FIG. 5: Positive energy vij1/2 states in 120S) under different box boundary conditions. 


The reason is that the centrifugal barrier for vfs/2 (J = 3) is much lower than that for vij3/2 (J = 6). The stability 
behavior of the state vfs/2 is presented in Fig. 3. The resonant energy and width are 0.870 MeV and 0.064 MeV, 
respectively. The wave function of this state is also well localized as shown in Fig. 4. 

Except the state vi,3/2, there is another neutron resonant state with | = 6, i.e., the vij,/2 state. As shown in Fig. 5, 
there is a less stable state lying at about 10 MeV. Although it shares the same centrifugal barrier with vi,3/2, vii1;2 
is much wider because its energy is much larger than that of vij3/2. In Fig. 6, one finds also that the wave function 
of Vi,3/2 oscillates very much even at r = 50 fm. The resonance parameters for this state are E, = 9.811 MeV and 
T = 1.275 MeV respectively. 

The energies and widths of single particle neutron resonant states obtained from the RSM calculations are sum- 
marized in Table II. We also calculate these resonances using the parameter set NL3 for the Lagrangian density in 
the RMF model and compare the present results with those from the ACCC approach and the scattering phase shift 
method [45] in Table II. In Fig. 7, the comparison is also made in a planar E-I plot. For the two low lying resonant 
states, vfs/2 and vi,3/2, the three methods give consistent results both for the energy and the width. Although similar 
energies are obtained from these three models for higher resonances, vij1/2 and vj15/2, clear differences occur among 
the widths from the ACCC approach and the scattering method and the results from the RSM lies in between. 
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FIG. 6: The radial wave function for the neutron vij1/2 resonant state in 1206n. 
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FIG. 7: Energies and widths of single neutron resonant states in '°°Sn from different methods. RMF-ACCC, RMF-RSM, and 
RME-S represents results from the analytical continuation in the coupling constant approach, the real stabilization method, 
and the scattering phase shift method in the framework of the RMF model with NL3 parameter set. 


TABLE II: Energies and widths of single neutron resonant states in '?°Sn from different methods. RMF-ACCC, RMF-RSM, 
and RMF-S represents results from the analytical continuation in the coupling constant approach, the real stabilization method, 


and the scattering phase shift method in the framework of the RMF model. In the RMF-RSM calculations, both PK1 and NL3 
parameter sets are used. All quantities are in MeV. 


RMF-RSM (PK1) RMF-RSM (NL3) RMF-ACCC (NL3) RMF-S (NL3) 


vlj Ey F 

vfs/2 0.688 0.032 
vi13/2 3.416 0.005 
vi11/2 10.01 1.42 


vjı5/2 |12.865 12.97 1.10 
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IV. SUMMARY 


In summary, the real stabilization method (RSM) has been developed within the framework of the relativistic mean 
field (RMF) model. With the self-consistent nuclear potentials provided by the RMF calculations with the parameter 
sets PK1 and NL3 for the Lagrangian density, the Dirac equation for the neutron is solved in the coordinate space 
under the box boundary condition. By investigating the stable behavior of the positive energy states against changes 
of the box size, the resonant states are singled out. The RMF-RSM is used to study single neutron resonant states in 
spherical nuclei. As examples, the energies, widths and wave functions of low-lying neutron resonant states in !?°Sn 
are obtained. Since a very large box size is used, even wider resonances can also be found. These results are compared 
with those from the scattering phase shift method and the analytic continuation in the coupling constant method and 
satisfactory agreements are found. 
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